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tool for the science exploitation of such a revolutionary mission, whose main goal 
is to trace back star directions from within our local curved space-time, therefore 
providing a three-dimensional map of our Galaxy. The results are useful for a 
thorough comparison and cross-checking validation of what already exists in the 
field of Relativistic Astrometry. Moreover, the analytical solutions presented here 
can be extended to model other measurements that require the same order of 
accuracy expected for Gaia. 
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To fully exploit the science of the Gaia mission (ESA, [I]), a relativistic astrometric 
model is needed able to cope with an accuracy of few /xas for observations within the 
Solar System. 

Gaia acts as a celestial compass, measuring arches among stars with the purpose 
to determine their position via the absolute parallax method. The main goal is 
to construct a three-dimensional map of the Milky Way and unravel its structure, 
dynamics, and evolutional history. This task is accomplished through a complete 
census, to a given brightness limit, of about one billion individual stellar objects. 

Since the satellite is positioned at Lagrangian point L2 of the Sun-Earth system, 
the measurements of Gaia are performed in a weak gravitational regime and the 
solution of Einstein equation, i.e the space-time metric, has the general form 

5a/3 — Vol/S T “b G (/l ^ , (1) 

where <C 1 and \dihai 3 \ ^ 1 can be treated as perturbations of a flat space-time 
and represent all the Solar System contributions to the gravitational field. Their 
explicit expression, however, can be described in different ways according to the 
physical situation we are considering. This means that, for the weak-field case, hafs 
can always be expanded in powers of a given smallness parameter e, as 

“b ^a/3 “b h(^jS -b * * * , 

(1) (2) (3) 

where the underscript n indicates the order of e. This expansion is usually made 
in power of the gravitational constant G (post-Minkowskian approach) or in power of 
1/c (post-Newtonian approach) both approaches coinciding inside the near zone of the 
Solar System [5] . While the post-Minkowskian formalism is better suitable outside the 
near-zone of the Solar System, the estimates performed inside this zone are sufficiently 
well supported by an approximation to the required order in (w/c) ~ e which amounts 
to about 10“^ for the typical velocities of our planet. Moreover, for the propagation of 
the light inside the Solar System, the sources of gravity should be considered together 
with their internal structure and geometrical shape. This is particularly true when 
the light passes close to the giant planets. In other circumstances it is an unnecessary 
complication to consider the planets different from point-like objects especially when 
the model is devoted to the reconstruction of the stellar positions in a global sense. 
At the microarcsecond level of accuracy, i.e. (u/c)^ ~ e^, the contribution to the 
metric coefficients of the motion and the internal structure of the giant planets should 
be taken into account, in particular if one wants to measure specihc light deflection 
effects, as for example, those due to the quadrupolar terms. 

The scope of this paper is to present an analytical solution for a null geodesic of 
the metric Q consistently with the requirements of Gaia’s astrometric mission and 
according to the RAMOD framework mm- RAMOD uses a 3-bl description of the 
space-time in order to measure physical effects along the proper time and in the rest- 
space of a set of fiducial observers according to the following measurement protocol 
0 : 

i) specify the phenomenon under investigation; 

ii) identify the covariant equations which describe the above effect; 

iii) identify the observer who makes the measurements; 

iv) chose a frame adapted to that observer allowing the space-time splitting into 
the observer’s space and time; 
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v) understand the locality properties of the measurement under consideration 
(namely whether it is local or non-local with respect to the background curvature); 

vi) identify the frame components of the quantities which are the observational 
targets; 

vii) hnd a physical interpretation of the above components following a suitable 
criterium; 

viii) verify the degree of the residual ambiguity, if any, in the interpretation of the 
measurements and decide the strategy to evaluate it (i.e. comparing to what already 
known). 

The main procedure of the RAMOD approach is to express the null geodesic in 
terms of the physical quantities which enter the process of measurement, in order to 
entangle the entire light trajectory with the background geometry at the required 
approximations. Then, the solution is adapted to the relevant lAU resolutions 
considered for Gaia [B]. 

Solving the astrometric problem turns out to compile an astrometric catalog with 
the accuracy of the measurement model. Indeed there exist several models conceived 
for the above task and formulated in different and independent ways iW and 
references therein). Their availability must not be considered as an “oversized toolbox” 
provided by the theoretical physicists. Quite the contrary, they are needed to put the 
future experimental results on solid grounds, especially if one needs to implement 
gravitational source velocities and retarded time effects. From the experimental point 
of view, in fact, modern space astrometry is going to cast our knowledge into a widely 
unknown territory. Such a huge push-forward will not only come from high-precision 
measurements, which call for a suitable relativistic modeling, but also in form of 
absolute results which need be validated by independent, ground-based observations. 
In this regard, it is of capital importance to have different, and cross-checked models 
which exploit different solutions to interpret these experimental data. 

For the reason above, inside the Consortium constituted for the Gaia data 
reduction (Gaia CU3, Core Processing, DPAC), two models have been developed: 
i) GREM (Gaia RElativistic Model, [7]) baselined for the Astrometric Global 
Iterative Solution for Gaia (AGIS), and ii) RAMOD (Relativistic Astrometric 
MODel) implemented in the Global Sphere Reconstruction (GSR) of the Astrometric 
Verification Unit at the Italian data center (DPCT, the only system, together with 
the DPC of Madrid, able to perform the calibration of positions, parallaxes and 
proper motions of the Gaia data). RAMOD was originated to satisfy the validation 
requirement and, indeed, the procedure developed can be conceived to all physical 
measurements which imply light propagation. 

Sectionj^lists the notation used in this paper; sectionj^is devoted to the definition 
of the mathematical environment needed to make the null geodesic explicit at the 
desired accuracy. In particular, in order to fully accomplish the precepts of the 
measurement protocol and to isolate the contributions from the derivative of the 
metric terms at the different retained orders, we introduce a suitable classification 
of the RAMOD equations. In section]^ we set the appropriate approximations which 
permit the analytical solution of the astrometric problem. In section we show the 
specific solution for the light deflection by spherical and non-spherical gravitational 
sources. In section hnally, we deduce the analytical solution of the trajectories of 
the light signal emitted by the stars and propagating through the gravitational field 
within the Solar System. In the last section we summarize the conclusions. 


Relativistic Astrometry for a Gaia-like observer 

2. Notations 


4 


• Greek indeces run from 0 to 3, whereas Latin indeces from 1 to 3; 

• ” : partial derivative with respect to the a coordinate; 

• scalar product with respect to the euclidean metric 

• ” X ”: cross product with respect to the euclidean metric bij ; 

• tilde-ed symbols ” refer to quantities related to the gravitational sources; 

• repeated indeces like Uav'^ for any four vectors means summation over 

their range of values; 

• ”[ajS\"'. antisymmetrization of the indeces a,/3 

• V: covariant derivative; 

• Vq.: a-component of the covariant derivative; 

• V(/): spatial gradient of a function / 

• we use geometrized unit, namely G = 1 and, c = 1. 

3. The RAMOD equations for Gaia 

The basic unknown of the RAMOD method is the space-like four-vector £“ , which 
is the projection of the tangent to the null geodesic into the rest-space of the local 
barycentric observer, namely the one locally at rest with respect to the barycenter 
of the Solar System. Physically, such a four-vector identifies the line of sight of the 
incoming photon relative to that observer. 

Once defined the equations of the null geodesic takes a form which we shall 
refer to as master-equations. Neglecting all the O(h^) terms, these read ITU]: 

_ 1 

— — tGhoj^i — -hoo,o = 0, (2) 

1 ___ 1 \ 

_ _ - hoo,r) + tP - 2 ] (3) 

P {^fc0,2 2^00, fc ^ Phoi,0 “1“ ^A:0,0 — 3- 

Here a is the affine parameter of the geodesic and 

d<j = dt + 0 {h ), (4) 

whereas t is the coordinate time. 

In order to solve for the master equations one should define appropriate metric 
coefficients. To the order of e^, which is what is required for the accuracy targeted 
for Gaia, one has to take into account the distance between the points on the photon 
trajectory and the barycenter of the a-th gravity source at the appropriate retarded 
time together with the dynamical contribution to the background metric by the relative 
motion of the gravitational sources. More specifically is the retarded distance 
defined as 

(5) 

where d is the parameter of the (a)-th source’s world line. The retarded position of 
the source is fixed by the intercept of its worldline with the past light cone at any 
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point on the photon trajectory. However, the retarded time t' = t — rt^a) and the 
retarded distance r(^g^ are intertwined in an implicit relation which would prevent us 
to solve the geodesic equations. Nonetheless, we show that it is possible to write an 
approximate form of the metric which retains the required order of accuracy of e^, but 
where the dependence from the retarded contribution is simplified. 

By using the Taylor expansion around any a' to the first order in e, we get for 
each source: 

a;*(cr) « + u*((t')(tt — it'), (6) 

which allows to rewrite the retarted distance as 


r\a) = x^a) - x^(a') « x"(a) - x^(a) + v"(a')(a' - a), 

i.e. 

r(a) =r"(cr,a)-hF(a,a')-i-0(v^), (7) 

where we set P((7, d') = x^(d') — x^(a'). Nevertheless (a' — a) is again proportional 
to the retarded distance as measured along the source world line. In fact, considering 
the tangent four-vector of the source world line 


where is the a-component of the spatial four-velocity of the source relative to the 
origin of the coordinate system and defined in the rest frame of the local barycentric 
observer u, the interval elapsed from the position of the source at the time and that 
at t is 


a' — a = ATu = —UaAx°‘ 

= — r]ai3U°‘Ax^ -j- 0{h) 

« Aa;° — dijV^Ax^ 0{h), (8) 

where Aa:“ = x°‘{a{t')) — x°'{a{t)). To the first order in v, we have along the generator 
of the light cone 


A.» = r(<,,i)+<!A 4 P^+ 0 (A 

r(cr,cr) 

then we get the following approximate expression for ([^: 

a)(a) « a) - P(CT')r(CT, d) 


(9) 

( 10 ) 


or 

r"{cr,^'){a) « r(cr,d)[n'(cr,d) - d'(d')], (11) 

where n*((T, d) = r*/r. This is equivalent, to first order in d, to the distance found in 
m and entering the expression of the metric, i.e. 

r{a,a')(^a)^r-r-v. ( 12 ) 

The choice for the perturbation term of the metric has to match the adopted 
retarded distance approximation and the fact that the lowest order of the h terms is 
and the present space astrometry accuracy does not exceed the level. 
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Then, for our purpose, a standard suitable solution of Einstein’s equations in 
terms of a retarded tensor potential mm, which can be further specialized as the 
Lienard-Wiechert potentials m is 


= E —+o 


hoi — — 


Ga) 

4.M(a) 

Ga) 


Pi{a) + O I 


(13) 


hij - + O (e'^) , 

a 'G) 

where Al(a) is the mass of the ath gravity source, r((j) is the position vector of the 
photon with respect to the source, (3^ = x^q = (1—/ioo/2)u*(d)+0(/i^) is the coordinate 
spatial velocity of the gravity source. 

Note that the time component of the tangent vector to the source’s worldline [3] 
is 


dt 


while that of the local barycentric observer is 
u° = ^ = l + ^+0{h^). 

do' 2i 


(14) 


(15) 


Then from (|14|) and (|15|) we derive the following relationships in the linear 

(16) 


approximation 


(id = dt(l-^-y) + 0(/i^) + O (e^) , 


and 


da = dt ^1 — + 0{hf) + O (e^) . (17) 

Within the approximation the perturbation of the metric transforms as 


doo = 2^^ i^')^(l + n(„) • V(„)) + C> (e'^) 

a 


, ^ -^{a) 

l^Oi 4 y ^ 


M 


r(a) 
(“), 


0{e^) 


(18) 


a '(“) 


or, by simplifing the notation 


^00 — — 2 A4(a)-(1 + n(a) • V(£j)) + 0(e‘^) 

a a 

hoi = - 2hvi ~ -2 /i(a)'Ci(a) + O(e^) 

a 

hij — hSij + ©(e^), 


(19) 


In the follows, ?'(„), unless explicitly expressed, will indicate a function with 
arguments a and d. Moreover, to ease notation we drop the index (a) wherever it 
is not necessary. 
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Let us consider a space-time splitting with respect to the congruence of fiducial 
observers u in the gravitational field of the Solar System [3]. The field equations 
can be rewritten in terms of the shear, expansion and vorticity of the congruence u 
[see [12] • For our purpose it is enough to consider only the expansion term and the 
vorticity [see fTO]. 

The master equations and (§ are obtained by retaining the vorticity term at 
least to the order of O{hoi), and the expansion to the order of O{dohoo) and O{dohoi). 
In the case of a vorticity and expansion-free geometry, the RAMOD master equations 
are named RAMOD3 master equation [3| 

-\- PP — -hij^kj + ^^'^Phoo.i — 2^oo,fc + O {hf') = 0, (20) 


where P = 0. Taking into account that PP6i 
(201 can be reduced respectively to: 

0{hf) 


dP 


1 , 


v){Php + -h^o 


1-1-0 (e^), equations ([^, and 


( 21 ) 


da 


+ h,k - \php - 2{P i)h,k + 2v'^{Php - ~PP{vP),o 


( 22 ) 


+ (?)'=h),o + 0(^2) (23) 

and for the static case {happ = 0 and z; = 0) 
c\P 3- 

-— + -P{fhp-hu + 0{h^)=Q. (24) 

dcr 2 

In order to fully accomplish the precepts of the measurement protocol, it would 
be useful, to isolate the contributions from the derivatives of the metric at the different 
retained orders. 

This allows us to classify the master equation as follows: 

• RAMOD3a (R3a), the spatial derivatives of the metric are considered while /ipi 
are neglected 

r\P 3 _ 

= --i\Ph,i) + + 0{h'^) = (R3a)'= (25) 

dcr 2 

• RAMOD3b (R3b), the spatial and time derivatives of the metric are considered 
while /iQi are neglected 

dP ^ 

— = -ho + 0{h^) = (R3b)° (26) 

da 2 

fifk 1 

— = (R3a)'= - -Ph 0 -f 0{h^) = (R3b)*^ (27) 

da 2 

• RAMODda (R4a), the spatial derivatives of the metric are considered including 
hoi 

dP 

= (R3b)° - 2{£ ■ vPPhp + 0{h^) = (R4a)° 


(28) 
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= (R3b)'= - 2{e ■v)hk + i) 

da 

+ 0{h^) = (R4a)'= (29) 

• RAM0D4b (R4b), the spatial and time derivatives of the metric are considered 
including hoi 

— = (R4a)° + 0{h^) = (R4b)° (30) 

da 

rjfk 

— = (R4a)'= - 2tt{vih) o + 2{v^h) o + 0{h^) = (R4b)'= (31) 

da ’ 

The implementation of RAMOD models and the need of testing them through 
a self-consistency check at different levels of accuracy, will benefit form this explicit 
classification. Beside this new classification of the RAMOD master equations, it is 
clear that the solutions call for an explicit expression of the metric terms. In general, 
for any integer m: 

(r)™ = mr*(r)’”-2, (r)™ =-to(/ 3 • r)(r)™-2. (32) 

From the last computation one could expect that in the case of mapped 
trajectories for RAMOD3-hke model [see IHD] the term a;*g should be retained, since 
each mapped spatial coordinate depends on the a value of the local one-parameter 
diffeormorphism. In this respect, note that the null geodesic crosses each slice S{t) at 
a point with coordinates a:*(cr(t)), but this point also belongs to the unique normal 
to the slice S{t) crossing it with a value a = a{x^,t) which runs differently for any 
spatial coordinate and therefore does not coincide with the proper time of the local 
barycentric observer. Therefore 

dax'^ = {dax\a)){dQa[x\X),t)) = 0 - 
Now, by using the retarted time approximation we get: 
ij-)~r{a,a') = - -b n^a) ■ + 2n(a) • v) - v^] + 0{h) (33) 

and 

(r);)j'"(cr,d') = mr(-pi(n(a) •/3)(l-bn(a) •/3)”"-b 0(/i). (34) 

Let us indicate the photon impact parameter with respect to the source position 
as: 

(35) 

and, for sake of convenience, let us denote also: 

(36) 

Finally, according to the previous derivatives, making them explicit, and denoting 
rd = rYr, the master equations assume the following expressions, valid up to the 
order: 

• RAMOD3a: 

^ = 2^ ^ ■ ^') - (1 + 2^(.) • i) + 

+0{v^) + 0{h^), (37) 

where in case of zero velocity we recover the static RAMOD recorded as follows 
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RAMOD3s(R3s): 

_ 9 -^(a) 

da ^ 


'(a) 




Similarly for the other classification items we have: 
• RAM0D3b: 

^=11 + 0{h^) 


(a) 


^ = (R3a)'" - X! 


da 

• RAM0D4a 
di° 


(a) 


M 


da 


= (R3b)° + 4^ ^ [{I- i){i- n(,))] + 0{v^) + 0(^2) 


'(a) 


M( 


^ = (R3b)'' +4^ —2^ ^ (^(a) X ^)]* + O(v^) + O(h^) 

CLCT V / \ 

a {a) 


(38) 

(39) 

(40) 

(41) 

(42) 


It is clear that, to the order of e^, we do not need to include the time derivative of 
liQk since these are at least of the order of e'* and should be neglected. Therefore, we 
do not consider the part of the RAMOD4 equations which contain the time derivative 
of hok and the second order velocity contributions. 


3.2. The case for an oblate body 

Now, let us consider the a-th source and define 

h(a) = ‘2M(^a)h(a) = 2Af(o); 


Ta) 


1 Jrr 


Go.) 


Pm(c0s6»(a)) 


(43) 


which means to take into account the mass multipole structure of the a-th body where 
Pm are the Legendre polynomials, Ma the mass of the body, Ra its equatorial radius, 
0(a) the co-latitude, and Jm the coefficients of the mass multipole moments. With 
this choice our considerations are confined to the case in which the object’ ellipsoid 
of inertia is an ellipsoid of revolution and the directions of the spatial coordinate axes 
coincide with those of the principal axes of inertia [13] . 

A rigorous treatment of a n-body multipolar expansion should take into account 
the different orientation of its axis of symmetry. However, this contribution decrease 
so quickly that at any accuracy currently attainable it turns out to be an unnecessary 
complication, since just one planet at a time would give a detectable effect. 

The derivatives of the metric coefficients, with retarded time approximation, have 
the following expressions: 


^{a),k I 2 4“ 271(a) ‘ H) T 2 



1 Jrr 


R. 


(a) 


r(a) 


Pm{cOS 0(a)) 


m(l -h n(a) • u)™Pm (cos 0(a)) 


Go) - 




- Pm(cOS0(a)).fc 


(44) 
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h(a),0 = { ^ I (1 + 


(o 


1 Jrr 


R(. 


Ga 


Pm(cOs6»(a)) 


1 


H-(1 + n(^a) ■ v) 

Ga) 


-Y.Jr, 


i?, 


(a) 


r(a) 


(P„^(cOs6»(a)),0 


+ mPn^(cOs6>(a)) ^^°2 ^ (1 + »(a) ’ '5)" 

^(a) 


(45) 


where n® = r®/r . A general n-body solution should include the multipolar structure 
of the sources. Nevertheless, according to the current astrometric accuracy and for an 
oblate body, the quadrupole approximation can be considered enough [see |7j . If we 
omit the higher multipole moments and restrict ourselves only to m = 2, denoting by 
s^ the axis of the sources which is normal to the source equatorial plane, Pm{cosd(^a)) 
is approximated as: 


P 2 (cos 61(a)) 


3(5(a) ■ r(g,d-'))^ _ 1 
2r(cr, ct')^ 2 


3 

2 


(^(a) ■ ^(a) 


- S(a) • ■5)^(1 +n(a) • V)'^ 


1 

2 ’ 


(46) 


which to first order in 0* becomes 

3 1 

P2(cOs6(a)) = -(S(a) • n(a))^(l + 2n(a) • v) - 3(S(a) • n(^a)){S(a) ' v) - ^ + 0{v‘^)] (47) 

therefore 

^(a),k 2 2^(a) ‘ ”6) “f ~~2 d” J^R 


3n 


(a) 


(a 


(a) 


'(a) 


5(s(a) •n(a))^ 


■ 9 ( 5 ( 0 ) ■ '^(a)) 


1 2 
(f^(a) • i') - 2 - 3(s(o) • ri(o))(s(o) • v) - -(n(o) • v) 

+ [1 - 3(S(a) • n(a)) ] + ^^[(S(o) • v) - (S(o) • ri(a))(l + 3n(a) • u)] 

^^(a) ^(o) 

+o(^"), 


(48) 


and 


^(a),0 — 


»T-(a) 


(o) 


3J2P^a) 

J.4 

^(o) 


»T-(a) 


■ U 5 , 


D)(. 


(o) 


1(a)) 


+ 0(0"). 


(49) 


By taking into account the target accuracy of Gaia [see [71 [14], the velocity 
contributions for an oblate body should be neglected. However, for sake of consistency 
and completeness with the assumptions adopted in this work, neglecting, a priori, 
terms which are part of the solution is not justified, even if the application to Gaia 
will surely dismiss many of them. Probably a Gaia-like mission that achieves few sub- 
microarsecond in accuracy will benefit of these analytical contributions, especially in 
regards of a cross-checking comparison between different approaches. In this case it 
would be better to consider a metric which properly contemplate all the complexities 
of a non-spherical gravitational body; that, at the moment, is out of the scope of the 
present paper and deserves a dedicated work (see, e.g., m)- 

Therefore, the RAMOD reduced master equations which take into account the 
quadrupole structure for the a-th single source finally become: 
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RAMODSaQ (R3aQ) 

2J2i?LX(a) 


da 


= (RSa)'' + 


(a) 


“(a) ^ 


5(s(a) •n(a))^ 


1 2 

+9(S(a) • n(a))^(n(a) ’ u) “ 2 “ 3(S(a) • n(^a)){S(a) ' v) “ 3 (^(a) ' v) 


dt- 




[1 - 3(s(a) • n(a))^] + 3 


\jk 


[0 


(a) ■ 


-(S(a) • ?T-(a))(l + 3(n(a) ’ u))] } + 0(v^) + O(h^) = (R3aQ)'=, 
where d* = sf^^ S(a)). 

In case of velocity equal to zero, RAMOD equations become: 




da 


(a 


'k ^He-na) 
a(a) ^ 


[ 5 (s(a) •^^(a))^- 1 ] 


+3 


Sg) 


(s(o) • n(o)) I + 0{h ) = (R3sQ) 


RAMOD3bQ (R3bQ) 


_K 

dcr 


= (R3b)° + 




2rf 

[a 


'^{a) ' ^ , ,2 

-^- ‘ V){s^a) 'n^a)) 


+O(0^) + O(/i^) = (R3bQ) 


dJU 

da 

RAMODdaQ (R4aQ) 

d£9 , 4M(a)J2Rj 

da 


= (R3aQ)'= - ^'=(R3bQ)° + 0(1)2) ^ (R3bQ)'‘ 


= (R3bQ)°- 4 - ■v)\-{i- Ua) [5(S(a) • n(a))2 - l] 

'^(a) I ^ 


-3(£ • Sa)(S(a) • «(o))} + 0(1)2) ^ (R4aQ)° 


(50) 


(51) 


(52) 

(53) 


(54) 


_K 

da 


= (R3bQ)'= - 


^(a) 


\{e X Ua X D)'= [5(S(a) • n(a))2 - l] 


+3(£ X D X Sa)''(S(a) ' ^(a))} + 0(1)2) ( 9 (/i 2 ) = (RdaQ)'". 


(55) 


4. Light propagation through the Solar System and parametrized 
trajectories 

When a photon approaches the weak gravitational field of the Solar System, heading 
to a Gaia-like observer, it will feel the gravitational field generated by the mass of 
the bodies of the system while it will be rather insensitive to the contribution to the 
field due to their own motion. If one compares the scale of the Solar System and the 
photon crossing time through it - approximately 10 hours in total- the gravitational 
field of the Solar System cannot significantly change in a dynamical sense during such 
time, to the point that the source velocity can be considered constant all along the 
photon trajectory. This last remark facilitates the solution of the RAMOD equations. 
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Let us make explicit the vorticity of the congruence u: 

= Pp“(u)P^(u)V[„u^] 

= ^ [pUa] + U[pUa] (56) 

where -P“(u) = Sf +u“itp is the metric induced on each hypersurfaces of simultaneity 
of u. Considering that = —1, and , we deduce: 

^0[p^fT]^00 “t“ ^[p^a]0 “t *^0 (^0[p^cr]o) 

which implies 


Woo = 0 

Woi = 0 

= d^ihj'^Q. 


W: 


(57) 


Taking into account the metric (19), equations (57) show that if we want a vanishing 


vorticity we have to choose hV x u + V(/i) x u = 0, which is satished if the velocity 
of the source is zero, i.e. a static case, or is constant, which corresponds to the case 
remarked above about the variations of the gravitational fields of the Solar System as 
regards the photon crossing time through it. 

Now, within the scale of a vorticity-free geometry, from the Frobenius theorem, 
the space-time can be foliated and one can always map the whole geodesic onto 
the hypersurface of simultaneity of the local barycentric observer at the time of 
observation. In this case the mapped trajectory can be expressed in a parametrized 
form with respect to the centre-of-mass (CM) of the gravitationally bounded system 

m- 


/ tdr, 


(58) 


where 


• C is the impact parameter with respect to the centre-of-mass of the gravitationally 

bounded system, i.e., = 0-|-(!l(/i) at the point of the closest approach with 

modulus I = 

• T = u — a, being a the value of geodesic parameter at the point of the closest 
approach. 


Furthermore, if we approximate the quantity P in terms of small perturbations 
with respect to the unperturbed light direction t'p': 

r = + (5r)2-b... ( 59 ) 

we note that the term {5t) can be neglected in the master equations and ^ 
because of the order of 0{h) (see also [H]). This implies that equation (58), for our 
purpose, can be approximated as: 


x* = + (4 + 5t)f + O (h^) , (60) 

where 6G is of the order of the deflection, i.e. e^. 

Finally, after these assumptions, the approximated retarded distance can be 
parametrized as 


r*(cr, a') « -I- (£^-l- 6t)T — 3A[d) — v^{d')r(a, d). 


(61) 
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and the distance can be reformulated as 

r*(cr, ct') « r*(f) — (62) 

where 

r*(f) = C + (4 + Se'^)r - (d) = fp + + «'=)f, (63) 

namely the relative distance on the slice at the time of observation without the 
contribution of the source velocity, while 


= C-i" 


(64) 


is the relative distance between the point of maximum approach of the photon to the 
CM and the source centre-of-mass. In figure Q it is sketched the relationship among 
these vectors on the slice corresponding to the time of observation. 

By using the scalar (•) and vectorial (x) products and the parametrization, the 
impact parameter ( |36[ ) with respect to the source becomes: 

d(f)'= = [£x (r(f) X £)]'=. (65) 


Note that it does not depends on r by definition, so it coincides with 

d’; = [i X {fp X e)]>^ (66) 

or, with the approximation ( [5^ , 

rfp=^p-4(^p-^», (67) 

with modulus dp = y/dp ■ dp. Moreover, we can assume that f/dp « 1/ tany, where 
X is the angle between the directions at the observer towards the point of maximum 
approach and the centre-of-mass of the source. In this case it is 


r(f) 



dl 


sin -I- 2 


■ Xp) 
dp 


cos X sin x- 


( 68 ) 


5. Light deflection by spherical and oblate spheroid gravitational sources 
with constant velocity 


Condition (57) constrains the solution of the geodesic equation to specific 
circumstances. In particular we do not consider equations Rda*^ and RdaQ’^ since they 
derive from the terms which are null because of our physical assumption that 

the sources move with constant velocity. As far as the light deflection is concerned, we 
expect that the velocity contributions become relevant in affecting light propagation 
in the case of close approach when general relativistic effects become of the order of 
Gaia’s expected accuracy together with the multipolar structure of the source. In this 
section we proceed to make the solutions for light deflection explicit according to the 
RAMOD classification. 
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Figure 1. Mapped trajectories and relative source positions 

on the slice at the time of observation (tq); is the impact parameter with 
respect to the centre-of-mass CM of the gravitationally bounded system, the 
unperturbed light direction, whereas i* that one at the star, and x the angle 
between the directions at the observer towards the point of maximum approach 
and the centre-of-mass of the source. 



5.1. Monopole contribution without velocity: R3s case 


This solution includes only the spatial derivative of the metric as in the static case. 
Then equation (38) can be integrated as follows: 








dr 




(69) 


namely: 






dr 1 ■ 


i.e. (see appendix A for the list of the solved integrals) 


a 

which reduces to 


1 


1 




rdr 


1 


0{h^ 


At 


R3s 


1 if 

1' 

1^2 

r 


dl 


Kn-rp)] 


0{h^). 


(70) 


(71) 


(72) 
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Formula (72) can be easily converted in the one found by Klioner (2003). If one 
considers a source at infinity = cr^): 


Ai'kss « 2 ^ M I - ^ -^Jl + a.n{n)]'j+0{h^). (73) 

Then, if the observable shift w.r.t. to the direction at infinity is (5 ct^ = c“^[crx Aix cr]^, 
we get 

Jcr'' « 2 ^ I - ^ [1 + cr • n{To)] 1 + O (li^) . (74) 


5.1.1. Monopole light deflection due to one body. Let us apply equation (38) to the 
computation of the light deflection due to one body, for example, the Sun, and check 
if the well known formula is recovered. 

With respect to the local barycentric observer u°‘ = e^Sg, the total light deflection 
is given by the modulus: 

P{u)apA.I.°‘ Afh = {gap + UaUp)A£°‘Ai^. (75) 

We have to consider only the euclidean metric 6ij, since A£^ is of the order of h, 
moreover it is already projected with respect to and in the static case P = 0; 
therefore, we limit the integration to the following expression 




1 


-dr — 


e 


1-00 [e + ^ 2 ) 3/2 


-df 


_2 7-oo (C2+f2)3/2 
which yields the well known solution: 

Al= 

We can also check the validity of our assumption by limiting the integration at the 
observer. Considering the source at infinity, from (68) 


,(76) 

(77) 


A^l3. = 2E-^Wl- 


4" 


smx 


1 + sinx2 + 


^ cos X sm X 


dj 


(£{/ ■ fp + To) sin X 


+ 1 


dpjl + sin x2 + 2^^^ cos x sin x 


0(/i2 


(78) 


and in the limit x << 1, since To/dp = cosx/ sinx, we obtain 


A£^ « 


(. -2- +0(.2). 


(79) 

(80) 


Then the modulus in the case of one body is again: 

- 4 

At= -.M. 

These results check the validity of RAMOD master equation when it is applied 
to one single body and make us conhdent to proceed with the solution for a n-body 
system. 
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This solution includes the perturbations depending on the spatial derivative of the h 
term and the retarted distance approximation. According to the assumption of section 
l^we consider also the velocity of the source constant, i.e. not dependent on f. For 
sake of simplicity, let us denote: 


jjk ^ '^p) 


'p 


w-») 

-^v — o’ 


(81) 

(82) 


In the case discussed here we limit ourself to the integration of the following 
terms: 

Ai’ksa « Ai'kss + 2 E \ r S + • v)] 


rdr 


r'^o x-2 


T‘‘df 


Di r s 


(83) 


In addition to the static case, the solution of the R3a master equation is: 


A^l3a = A£l3,+2E 


M 


(a) 


dj 




-Dp{fp-v) 


■ n) 
r 


ok 


^vdp “ ■ ^) + ^ ^ ■ (^P ^ ^9')) 


arctan 


'O'- 


Dl{iri)-^{dp-v) 



■(fp • n)' 



r 



+ 0{h) + 0{v'^) 


(84) 


If the source is at infinity and in the limit of x << !> i-®- the observer and for 
grazing rays, these contributions vanish (see appendix A). 


5.3. Monopole and velocities contribution: R3b case 

This solution includes the perturbations depending on the time derivative of the h 
term according to the retarted distance approximation. Assuming the source velocity 
constant, we have to integrate: 


A^3{, 


dr 


fdr 






The solution is straightforwardly (appendix A): 

M^a) 

dp 


(85) 

( 86 ) 


A^Sb « E {(^P ■ ■ ^p)]? } 


A/'fe ~ —P^aP 

^'^R3b ~ ^^R3a ^p^’-R3b- 


( 88 ) 
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For the source at infinity, in the limit of the small angle x between the directions at 
the observer towards the source center of mass and the point of closest approach, from 
appendix A these contributions result: 


^^R3b " 

a 

to 

^^R3b " 

a 

K 

1^^ 


(£y ■ Up) 
2 


Ik ifp ■ v) 


(89) 

(90) 


5.4- Monopole and velocities contribution: R4a case 

The R4a^ equations include the perturbations depending on the spatial derivative of 
the IiQi term, i.e. those depending on the mass-current contribution. According to the 
constrain on the vorticity imposed by the physical assumption as discussed in section 
Q, we have to consider only equation R4a°: 



" A^gj, 


" ^^R3b 

i.e. 


A44a ^ 

" A^g^ 

A^i)4a ^ 

^ A£lg, 




rdr 

^3 


~ ^^R3b + ■ ^) {■ ^p) [(^9' • n)] - [{n ■ fp)]l° I 


(91) 

(92) 

(93) 

(94) 


Again for a source placed at infinity and in the limit of a grazing light ray, we get 
at the observer: 


A£?e4a ^ 

a 

A44a ^ 



dl 


dl 


{{fp-v) + A{egfv){egf-fp)} 


1-f 




Ik (rp ■ v) 

V 2 


(95) 

(96) 


5.5. Light deflection by a static oblate body: RSsQ case 


For sake of convenience, in this we drop the suffix (a). Considering that the quadrupole 
structure of the source may be relevant in affecting light propagation only during close 
approach, here and after, to compute the quadrupole contribution, we refer only to a 
single body. Moreover, as in section (5.2), we denote: 

(97) 


where s^ is the axis of the source normal to its equatorial plane. Equation (51) in 
function of the parametrized trajectory can be integrated as follows 


At 


dr 


dr 


R3sQ 


= 3J2R^M <-[Dl + 2(s • fp)D'^] / ^ + 5(s • rp^D^ / ^ 


(98) 
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+ 





+ 5(s • fp) [2(s • 1^)0’; 



+ Hs-i^)[is-i^)D^p 



54(s • 4)2 




+ o{h^), 


Table 1. List of the recurrent integrals 


Notation 

Corresponding expressions 

r(f) 

^f2 + 2(£p--rp)f+ f2 

C/0 

arctan[(^Q'- r)/dp\ 

h/0 

Log(r ■£(>•+ r) 

C/i 

Ijr 

C/2 

l/r^ 

C/3 

l/r3 

C/4, 

l/r4 

C/5 


C/5 

1/r® 

C/4 

(C'CC 

C/2 

(C' r)/r'^ 

C/3 

iC'C/C 

C/4 

(C'CIC 

C/5 

{C-r)lr^ 

C/0 

iC' r)/r® 

C/4 

(r ■ fp)/r 

C/2 

(r ■ rp)/r^ 

C/3 

(r ■ fp)/r^ 

C/4 

(r ■ fp)/r* 

C/5 

(r ■ rp)/r^ 

C/0 

(r ■ fp)/r^ 


namely (for the details of the solution see tableand Appendix A), 

= 6J2i?^ ^ {[/0/3 ]Cf + [f/0/3] r*’ + [f/0/5 ]+ [h/i] 1“ cl 

+ [^ 1 / 3 ]/ C^ + [-fi/ 5 ]/ Cq + [d2/3]/ + Cy [^ 2 / 5 ]/ C's } + 0{lr^), (99) 

If we consider a grazing light ray emitted by a star located at infinity, above equations 
reduce to: 

= 6J2i?"^ {2C|} + 0(/i2). (100) 

Similarly to the spherical case, let us compute the total light deflection close to 
one gravity source. In this case, we expect that our formula reduces to the available 
expressions known in the literature [laillilTellli]. For this scope we can assume: 

• dp = fp = the radial vector; 

• dp/dp = —n^ the unit radial vector; 

• 4 = the tangential vector to the line-of-sight; 

• V = dp the impact parameter, 

• and = {£ x n)^ the orthoradial direction. 
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Then: 


,M 


^^msQ = + 2(s • n)[s'= - t'^{s • i)]} + O{h‘^)(,101) 

and, by expressing the vector s^ as linear combination in terms of the orthonormal 
basis (t,n, m) 

= (s • + (s • ri)n^ + (s • m)m^, (102) 

we obtain the same formula deduced in Crosta and Mignard m- 

Ae^ = Ai>k,s + ^i%3sQ 


4M 

dn 


1 (1 - (s ■ - 2{s ■ nf) 


v!^ + 


(Pp 


{s ■ m){s ■ n)m^ > (103) 


5.6. Quadrupole contribution with velocity: R3aQ case 


Equation (50) in function of the parametrized trajectory becomes: 


^^R3aQ — ^^R3a + ^^‘r3sQ 


+ 2J2R^M 

3 


D^^+SDfis-v) 


dr 


- [-Dp(9(s-rp)(s-u) + 2(fp-u)) 


dr 


dr 


(104) 


-2^vis ■ fpf + 3£>^(s • fp)(fp • 5)] J —+ 27D^{s ■ fpf(fp ■i)J — 

- [Dp (9(s • l' 9 ')(s • v) + 2(£9'- w)) - i^{%s ■ fp){s ■ v) + 2(fp • v)) + 3ZlJ(s • i^){s ■ fp) 

r^o -j- 

{[fp ■ 7 j){s ■ l^) + {s ■ fp){itji ■ v))] / +27[s-fp)[Dl{[s-fp)[£^-i) 


iy\j 

fdf 

J.8 


+ 2(s • £g)[fp ■ i)) - 4(s • fp)[fp ■ W)] I ^ + [9^(5 • igjis ■ v) + 2&y[lyf)) 

1 


-iD^Xk- ((^ • 4) + (s • X) + 


+ 27 [{s ■ XD. 


+ (2(s • fp)(£ 9 '- w) + (s • Xi^p ■ v)) - £X ■ X + ((s ■ fp)X- v) 


+2(s-V)(rp-u))] 
+ (s-£9')(fp-u))] I 


a2 


T^dr 

+ 27[s ■ X [dX ■ mk ■ v) - 4 (2(s • rp)(4 ■ D) 

3 27 [£> • iXX-v)] r + OX) + OX). 


yjO 

f^df 


Namely, from table 1 and appendix A, we denote 

AIXq = A^aa + ^imsQ + {[/o/o];° X + [h,X; X + [rh/AY; X 

+ Yo/e] f O12 + [''■l^o/e] f + [71/2] f C'w + [A/4] ^ C'ls + [A/e] ^ C'le 

+ [72/4] CY + [72/6];° of,} + OX) + OX)- (105) 

For a source at infinity and in the limit of % << 1 the above contributions vanish as 
in the case of mass monopole moving with constant velocities. 
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From equations (53) we have to solve 


^^mbQ — + “iAi'hR^ 




(106) 


dr 


j —-5(s-fp) 


1 


(fp-5)(s-£p') + -(^p^-5)(s-fp) 


-5(s • iy) 


7i{s■^v){fp■v) + {^^■v){s■fp) 


f^dr 



[is ■ J + Oiv"^) + 0{h^) 

A£l 3 ,Q = Ai'k.^Q - £^(R3bQ)0 + Oiv^) + Oih^). (107) 

Expressing the above formulae in terms of the solutions of the integrals listed in table 
1 , these quantities read 


^^R 3 bQ - ^^R 3 b -I-^— {[4/3] f C'lg + ["^4/3] / C'go + ['^4/5] / C21 + [h/i] / 

C 22 + [ 11 / 3 ];“ C 23 + [/l/5];° C 24 + [/2/3]I“ C^25 + [/2/5]^“ ^26} + 0{v^) + O(/l2Xl08) 

^^R3bQ = ^^R3aQ ~ ^^‘^^R3bQ + 0{v^) + 0{h‘^). (109) 

For a stellar source at infinity and in the limit of light ray grazing the gravity 
source, we get: 


^^R3bQ — ^^R3b C 22 + 0{v‘^) 0{h?) (HO) 

'^^mbQ = ^^R3aQ ~ ^I^^R3bQ + 0(5^) + 0{h?'). (HI) 

With the same notations of the single static source, the above equations collapse to 

^^R 3 bQ = ‘^^R 3 b + {-(^ H)[l - 4(s • nf - (s • <)2] + 2(s ■ t)(s ■ n)(t ■ u)} 

+0(v^) + 0(h^) (112) 

^^R 3 bQ = ^^R 3 aQ + “^^ 3 ^ {i^i^ H)[l - 4(s • nf - (s ■ tf] - 2t'=(s • t){s ■ n){t ■ u)} 

+0(u2) + 0(h2) (113) 


5.8. Quadrupole contribution with velocity: R^aQ solution 


The integral expressions of equations (54) and (55) with the physical assumption of a 
source moving with constant velocity are: 


1 

2 


• Tp) -\- 3{s • £(/){s • Tp) 



dr 


^^R4aQ — ^^R3bQ ~ 4A1 (a) ^2 -R^q) (l^Q'’ 5) 
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■ 3(s • 


rdr 


rdf 




dr 5, . , r„/ TT N 

-7 + [2(s-^9') 


i^vrp) + is-fp)] / ^ +-(s-fp^) [(£ 9 '-fp)(s-fp')+ 2(s-fp) 


r'^’o ^2 


T^df 




T^dr 


+ 0{v^)+0{h^) 


^^R3bQ — ^^RiaQ + 0{v^) + 0{h^) 

Considering table 1 and appendix A, the solution reads: 


2M(a)J2R\ 


^^RiaQ — ‘^^R3bQ 


dl 


(114) 

(115) 


-^0/3. C '2 


C 23 + 


27 -I- ['^^/a] 

C29 + [A/i]/ C'so + [A/3] / Aai + [A/s]/ A32 + [A/a] / A33 + [A/s] / A34I 

+0(f2) + 0(f/ 


(116) 


^^fi4aQ — ^^R3bQ + A(f^) + 0(/l^) (117) 

For a source located at infinity and in the limit of the small angle x between the 
directions at the observer towards the body center-of-mass and the point of closest 
approach of the light trajectory, we obtain : 


^^R4aQ — ^^R3bQ 


4A4(„)J2i?/ 

dl 


-(A-^)A30 + O(7") + O(d2) 


(118) 


^^R4aQ — ^^R4aQ + 0(v^) + 0{h^). (119) 

Again, with the same conventions of the static source, for the time component we 
deduce finally for the light deflection: 

16M(a)Al?(„) 


‘^^R4aQ — ^^R3bQ + 


dl 


-{t • f;)(5 • n){s • t) + 0{v"^) + 0{h^). 


( 120 ) 


6. Photon trajectory 

The first integration of RAMOD equations gives the estimates of the deflection 
effects. The second integration solves the ray tracing for the photon emitted by a 
star and intercepted at the observer’s location, thus obtaining the quantity Ax“ = 
a;“(fo) — a;“(f*) at the different approximation levels. In general since the first member 
is (ix“/df — the integral to be solved in the case of monopoles with constant 
velocities became: 


= ^^Af + 2 ^ M^a) r ^o/idr - J // h/idr \ + O (h^) , (121) 


jfc /.f„ 


d2 

“p J T, 


^^R3a — ^^R3s 


^ dl 


-c^pvAd-w) / ^Io/2dT-I A/2df 
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D 


Dydp “ ^ ^ ^p) ■ (^P ^ ^9)) 


lo/odf 


pk 

^9'/ 


D;ierv)-^{dp-v) 


Wrj +0(z;2) +0(/i2) 


( 122 ) 


I (fp -0)^ J / 2 /ldf 

+0 (z; 2 ) + O (/i 2 ) 

^^R3b = ^^R3a ~ ^^^^R3b + O [v^) + O (h?‘) , 
and 


M 


^a;?j4„ = Aa;?j3j, + 4^-^(^p.-0)|(£p.-fp)y^ Ii/idr-J^ / 2 /idf 


+0 (v^) + O {h^) , 


^^R4a — ^^R3bi 


(123) 

(124) 


(125) 

(126) 

where If = £°‘(to) is the line-of-sight direction at the observer used also to solve the 
boundary value problem [18], and the relationship At = (£q'- Ar) holds. 

Then, the trajectory is composed by the following terms (see appendix A): 

Aa;^ 3 s =^o^'^ + 2^M(„) [Log(r-^p- + r)]!° - J [r]ll\ + O {h^) , (127) 


^^R3a — ^XR3 s + 2 ^ 


M( 


D 


dl 


dp{£(^-v)- 


arctan 


dvdp - M ^ ^ 


{£{^ ■ r) arctan ( 


'P / Jr* 
r ^ 


^k 

dldl + Y ((n X fp) • (rp x £ 9 -)) + ^(dp • i)(rp ■ £^) 


[Log(r)];; 


+0{h^)+0{v^), 


-I r* 


(128) 


Aa:?j 3 {, =l)(Af + ^M(a) |-( 4 -^) [Log(r • £ 9 -+ r)][r]T°| 

+0 (^ 2 ) + 0 (?; 2 ) 

^^R3b = ^^R3a ~ ^^^^R3b + O (/l^) + O(u^), 

and, 

^X°p.ia = - 4 XI -^(a) (V ■ [Log(i' • h + !■)] 

a 

^^R4a = ^^R3b + O ih?) + 0{v^). 


(129) 

(130) 

(131) 

(132) 
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Similarly, the inclusion of the quadrupole terms, at the first order in h and v, 
gives additional contributions to the following trajectories (see table 1 and appendix 
A for the detailed expressions of each term): 


Aa:« 3 ,Q 6J2R I ^2 


3d4 


Cl. Cl 




dl 


Ml 


[^2/3]/. + C”! Hf° - C'l - f - [^0/3] + C? ['^^0/1]/. + [tIo/s]^, \ ( 133 ) 


^^R3aQ — ^^R3a + ^^R3sQ + 


.M 


„fc _ A„fc , A^/c , , ) I 


dl 


k , Cl, Cl,{i^-rp) , 3 CI 2 


dl 


2dl 


8dl 


3Cl,{£g/-fp) Cl-J 3C|8 1 r -|T„ nk\ 2 ^j I (fik 1 f-ik\ 

-8d|-+ ^ + ^ ) + C, [r/o/oj, + (C,4 - dpC,) 


'P 

[Log(r)]r: + 1 {cl, + cl, + a, ■ rp)Clr) [I2/2] ^ ^ (c|o + 


4ci2 


4d2 




• ^p)C'i5 + + 


4d2 


4d2 


iCl3 + {i^-rp)Cl, + Cl,) [ly^l: 


(134) 


Aa^flshQ — Ax^jst, 


■ 2 AI f Ci9 2C2i(£p'- Tp) ^ 2C26\ tt I'j'o C'zi 

Ml 




3 . , [A/J?. - ^ [V.]?- 


(hP" -PhP —'^ 22 H 1 ° - C23 [ 4 /i]/.- Y [-^0/3]/* 




+ C 25 [f/o/r];:+^[f/v3j,. 

^^R3bQ = Ax|j3f, — 


(135) 

(136) 


and, finally, 

C _ A.„0 


,2M 


^^''R4aQ — ^^R3bQ 2 J2-R^ ^2 (^9'' 
1^28 r. ] To / C 29 \ 

+r^33 [^^/i].: + ^ [^^ 0 / 3 ]/. 1. 

^^R4aQ = £^^R3bQ- 


C 27 2C29(1'p'-T p) I 2 CI 4 ^ 

Hi 3d4 +^) 

C: 


[h/i\A, - ( [l'2/3]/. +C3o[r]l°. +C 31 [/o/i]/.-^ [Vs]/* 


(137) 

(138) 


7. Conclusions 


The Relativistic Astrometric MODel (RAMOD) is a mathematical tool conceived 
to model the astrometric measurements made by an observer in space. Since its 
original purpose was to address this problem for the ESA Gaia mission, whose final 
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astrometric accuracy requires the physical model to be accurate at microarcsecond 
level, RAMOD had to take into account the general relativistic corrections due to the 
bodies of the Solar System. Despite the apparent straightforwardness of the task and 
the linearity of the metric given the weak gravitational field regime inside the Solar 
System, the solution of the inverse ray tracing problem, which allows us to reach the 
aim above, is rather intractable unless treated numerically, particularly if retarded 
time contributions need to be accounted for [4]. As far as RAMOD is concerned, the 
reason lies mainly on the fact that the main unknown of the differential equations is 
the observed direction as projected on the rest space of the local barycentric observer 
and represents locally what the observer measures of the incoming photons in his/her 
gravitational environment. This aspect transforms the geodesic equation into a set of 
nonlinear coupled differential equations which comprises also 

for the time component. The original version of the RAMOD model was 
therefore numerical and although successful in its applications with the inclusion of 
the relativistic satellite attitude [18) . it was hard to control and compare with similar 
astrometric models even with a comprehensive error budget for stellar positions m- 
Here we present a fully analytical solution of a system of differential equations 
up to the level everywhere in the Solar System, which is therefore able to assure 
a microarcsecond-level accuracy consistent with the precepts of the measurement 
protocol in General Relativity, and that can also be utilized under observing conditions 
more demanding than those of the Gaia mission. 

The analytical solution is general enough to be applicable to other missions 
conceived to exploit photon trajectories and extends the analysis of the trajectory 
perturbations due to gravitating sources with a non negligible quadrupole structure. 
While the retarded time approximation adopted here and the solution for the static 
cases recover the results obtained by similar astrometric models (as proved also in 
[To]'), the solutions including the constant velocity of the source give rise to different 
expressions that deserve to be carefully evaluated in a separated work as were done, for 
example, in [5], and |20j . At first glance, in fact, the presence of the time component 
for the local-line-of-sight has not been contemplated in other models and the 
components do not show complete coincidence. A proper comparison, both analytical 
and numerical, will precise the physical significance of this unexpected discrepancy as 
one carries on with the implementation process of the astrometric observables from 
which the relativistic astrometric parameters with their appropriate variance and, 
possibly, covariance values are deduced. 

Acknowledgment s 

The authors wish to thank Dr. Donato Bini and Dr. Andrea Geralico for the fruitful 
discussions and useful comments. This work was supported by the ASI contract 
1/058/10/0 &/1. 

References 

[1] Turon C, O’Flaherty K S and Perryman MAG (eds) 2005 The Three-Dimensional 
Universe with Gaia 

[2] Hawking S W and Israel W (eds) 1987 Three hundred years of gravitation 
(Gambridge University Press) 


REFERENCES 


25 


[3] de Felice F, Crosta M T, Vecchiato A, Lattanzi M G and Bucciarelli B 2004 
Astrophys. J. 607 580-595 

[4] de Felice F, Vecchiato A, Crosta M T, Bucciarelli B and Lattanzi M G 2006 
Astrophys. J. 653 1552-1565 {Preprint arXiv: astro-ph/0609073) 

[5] de Felice F and Bini D 2010 Classical Measurements in Curved Spaee-Times 
(Cambridge University Press) 

[6] Crosta M and Vecchiato A 2010 Astron. Astrophys. 509 A37 

[7] Klioner S A 2003 Astron. J. 125 1580-1597 

[8] Teyssandier P 2012 Classieal and Quantum Gravity 29 245010 {Preprint 1206. 
6309) 

[9] Hees A, Bertone S and Le Poncin-Lafitte C Phys.Rev.D 

[10] Crosta M 2011 Classical and Quantum Gravity 28 235013 {Preprint 1012.5226) 

[11] Kopeikin S M and Schafer G 1999 Phys. Rev. D 60 124002 

[12] de Felice F and Clarke C J S 1990 Relativity on curved manifolds (Cambridge 
University Press) 

[13] Bini D, Crosta M, de Felice F, Geralico A and Vecchiato A 2013 Classical and 
Quantum Gravity 30 045009 

[14] Le Poncin-Lafitte C and Teyssandier P 2008 Phys. Rev. D 77 044029 {Preprint 
0711.4292) 

[15] de Felice F and Preti G 2006 Class. Quantum Grav. 23 5467-5476 

[16] Epstein R and Shapiro I I 1980 Phys. Rev. D 22(12) 2947-2949 URL http: 
//link.aps.org/doi/10.1103/PhysRevD.22.2947 

[17] Crosta M T and Mignard F 2006 Class. Quantum Grav. 23 4853-4871 {Preprint 
arXiv:astro-ph/0512359) 

[18] Bini D, Crosta M T and de Felice F 2003 Class. Quantum Grav. 20 4695-4706 

[19] de Felice F and Preti G 2008 Classical and Quantum Gravity 25 165015 

[20] Bertone S, Minazzoli O, Crosta M, Le Poncin-Lafitte C, Vecchiato A and Angonin 
M C 2014 Classical and Quantum Gravity 31 015021 {Preprint 1306.2367) 

Appendix A. Expressions of the coefficients in the quadrupole deflection 

formulae 







REFERENCES 


26 


= - 1+4 


{s-fpf 


[-5(s • fp) - 3{fp ■ ds)] + fl 


" dj 


+ ^ {-(^'p • ^(/) + rp) ^ ^ [-(s • Tp) 


+3{fp ■ £(/){s ■ £i^) 


ifp-e<^fis-i>^y\ 2 D 


3 dl 


ci = d;\-- + 


1 2{s-rpY {f-p ■ £(/){s ■ £(;) [4 


^(s-f-p) + 


{rp-£^){s-£M (s-£t/)^r 


+ ^ |- 2 (^p-V - (s-^p) 

-D^(s-rp); 

Ce^ = ^(^-rpr, 


2^2 (s • ^9')(s • rp) (rp ■ £i/)^(s ■ £y) 


(-2(s • Tp) - (fp • ds)) 


Cl) = -f + is- k)D^s-. 


Cs = —^(s • fp)^(rp • ^9-) - y {(s • rp){rp ■ £^) + (s ■ lyf [d^ - (fp ■ £(/)'^] } 


Cg = yy |-[9(s • fp)(s • v) + 2(fp ■r)] + ^(s- fp) (fp • dy) + (fp ■ f9')[9(s • £f;)(s • v) 
+2(4- f)] - ^(fp • £i^)(s ■ £^)(fp ■ v)(s ■ fp) + ^(s • 4)[4(^p • ^9-)^ + fl][2(s ■ fp)(£^-v) 
+ (s • £(/)(fp ■ f)] + ^[2(’'p • + 3^p](s • ^9')^(^P • 4')(^9'-1')]| 


+ ^ (--)- 


• fp) + [3(s • fp)(fp ■ £^) - (s ■ £i^)[2(fp ■ £yf- + f\ 


+ 7:1 \ “(^p ■ ^P')[9(s ■ ^p)(« ■ r) + 2(fp • f)] + ^(Pp • £y)(s • fp) (fp • v) + 


[ 2 (fp-^»^+f^] 


— — 9f4f T • £n'')^ -h 7’^] — — 

[9(s-v)(s-£(^) + 2(£(^-i)] -+yp-+^[(^9'-1') ((s • ^p)^ - ^p(s • ^9')^) +2(s-fp) 


, 9 [ 2 (fp-fp^) 2 + 3 f 2 ] 

(s • V)(rp • v)]-- 


(fp • ■^ 9 ')(s • ■^ 9 ')[ 2 (s • fp)(£i^-v) + (s ■ £^)(fp • f)] 


Dl J 9 (s-fpf 9 (s ■ £(^)(s ■ fp)(fp ■ £t^) (s-f9')2[2(fp-f9')2+f2] I 

2d„ 1 Adi 2dl Adi I’ 


Cfo = - 27 i 4 p" (s-f» 2 (f>.f) 


(rp • ■^ 9 ')^ „ -2 


(fp-4)" , 2f2 + (fp.f»2l'l 274 


{(s-^v) 


+ 2fp + (fp • £^) [2(s • fp)(£,^-v) + (s • £(^)(fp ■ f)] + (s • £ 0 ) (£iyv) 


K3 



REFERENCES 


27 


17(fp • ^^)f^ - 28f^(fp ■ (.gY + 16(^p • 
3d2 




{fp-^g)fl\ ^ idlD^ 


(f 


H 


+ {(-V «)(s • V) - • kY 


+Y^{dl^{s • ^9')(s • v) + {^g■v)\ + [rl{fp • ■ v){s • fp) + (s • 


(fp • u)] + 9(s • u)[7fp - 20r2(fp • £9^)2 - 8(fp • } ; 


’p “-'pV 

<^12 = QD^dlf^is ■ lyY{iy • u) + {(s • 4)^(^9' • ^Ylirp ■ hf “ dlr^is ■ iy) 

[2(s • rp)(^ 9 '• w) + (s • l^)[rp •«)]}; 

Cia = 9Dpdl{s • 4')[2(s • ^p)(^p'- ^') + (s • ^p')(^p • v) + 3(rp • £p')(s • 

+9^ {(s • 4)^(4- ^')[^P - 5rfp(rp • 4)] - 3rfp(s • ^9')(?'p • 4)[2 (s • rp){I^-v) 

+ (s- V(rp • {;)] - (s • fp)dl[{s ■ fp){it/-v) + 2{rp ■ i;)(s •4')]} ; 


rk 


45 

9(s • fp)(s • v) + 2{fp ■ v) + ^{s ■ rpYirp ■ dy) + {fp ■ £p')[9(s • ■ v) 


+2{e(/-v)] -^(rp- V)(s • ^9')(rp • v){s ' fp) + ^ [4(fp • - {s ■ ^^){rp ■ £>^) 

({fp ■ W - [2(s • rp)(4 • v) + (s • 4)(^P • «)]I + |(s • w) + 


dl 


(V • v) 


[(fp • i))(s • ^p-) + (s • rp){£^-v)\ -|^[2(^P ’ ^dY + ^p](s ' ^Y)\ + {-(^p ’ V) 


34 


Ml 


45 - - 3 - 

[9(s • fp)(s • v) + 2{fp • {))] + ^(^P • ^P')(s • fpYifp • v) + [2(fp • £^Y + ' v){s ■ (f/) 

1 ( 4 - «)] - ^ ^ ^ ^ [( 4 - ((s • fpY - ^p(s • ^ 9 ')^) + 2 (s • fp){s ■ l^){fp ■ u)] 


2dl 


9(fp ■ ^y)(s ■ £y)[2(fp ■ + iPp] 

2dl 


[2{eg,-d){s ■ fp) + {s ■ £^){rp ■ u)] 


_^_^il_Y-rpY , 9(s • £9 ')(s • fp)(fp • ^p-) ^ (s-£p')2[2(fp-£p')2+r2] 


dl 


2dl 


Adi 


C'ls = ^ ■ ^pYs ■ v) + ^(rp • i) + ^(s • fp) [(s • fp)(fp • v) - {fp ■ £g) {{£g ■ v) 

£fy')^ -U 7*^ _ _ 

-^[2(s ■ fp)(£r f) + (s ■ IgWp • ^)] 


(s • fp) + 2(s • fp')(fp • v))] + 3(s • fp-)- ^ 


3i:»J, _2 3L»,\ ,, 94 

-3(s • fp)[4(fp • lp')2 + fl][{ig- f)(s • fp) + 2(s • Ig){fp • f)]} ; 



REFERENCES 


28 


Cf6 = 9i?p"(s-fp)2(fp-0); 

<^17 = ^ {[9(s • 4)(s • + 2(4- ^)]} + y {(^p • 4')[9(s • ^9')(s • + 2(^9'- w)] 

3£)fc 

-[9(s • fp)(s • w) + 2(fp • u)]} + ^[(s • ^9)iTp ■ v) + {{s ■ fp) 

- 

-(s ■ £»(fp. y))] + ■ i^rifp ■ ly)- 

C% = -9Dp {3(s • ■rp)[(s • fp){i(^-v) + 2(s • l'p')(fp • v)] - (s • i>/){fp ■ ^9')[2(s • fp)(i(^-v) 
+ (s • ■ w)] - 3(s • i^f{i^-v){fp ■ l^f] - 9£^{3(s • fpf{rp ■ v) + {s ■ rp){fp ■ I^) 

[(s • fp){i(/-v) + 2(s • i(^){fp ■ w)] + (s • e(/)ifp ■ ^ 9 ')^[ 2 (s • rp){it^-v) + (s • / 9 ')(rp • v) 

-(ip(s-£p')2(fp-4)(4-^)]}; 

(^19 = -^ {[2rfp - SPpdp + 3(fp • 4)^)](s • 4)^(4- ^')} ; 

C20 = -^(s • 4)^(4- ^')(^p • 4); 

C 21 = 3dp(s • 4') ^{rp ■v)is- £(/) + ■ v){fp ■ ds) ; 




REFERENCES 


<^29 = dl{s • ^p')[(rp • ■ £(^) + 2{s ■ fp)]; 

^30 = ■ dp) - • ^9') [2(s • rp) ((s • ii^)ifp ■ it;/) + 2(s • fp)) 

-{s-ig)^[fl + {rp-ltif]] ; 


C'ai = -(^p • ■^p') - 6(rp • s)(s • (-t/) + ' '®)(^p ' ^9')[2(5 ’ ^9')(^p • ^9') + 3(rp ' 

, ?^p + 3(rp • p;,oC' M _ (!'p•4)^(s•4)^ 

(s • • ^9') (5 • £p') + 2(rp • s)] 


3fi2 


3(i2 


(^32 = (s-fp)^(fp-£p'); 


C'aa = ^[1 + 6(s • £{/f]. 



